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Abstract 

We present a very quick and powerful method for the calculation of heat-kernel 
coefficients. It makes use of rather common ideas, as integral representations of 
the spectral sum, Mellin transforms, non-trivial commutation of series and integrals 
and skilful analytic continuation of zeta functions on the complex plane. We apply 
our method to the case of the heat-kernel expansion of the Laplace operator on a 
D-dimensional ball with either Dirichlet, Neumann or, in general, Robin boundary 
conditions. The final formulas are quite simple. Using this case as an example, we 
illustrate in detail our scheme — which serves for the calculation of an (in principle) 
arbitrary number of heat-kernel coefficients in any situation when the basis functions 
are known. We provide a complete list of new results for the coefficients B3, Bio, 
corresponding to the D-dimensional ball with all the mentioned boundary conditions 
and D = 3,4,5. 
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1 Introduction 



An important issue for more than twenty years now has been to obtain explicitly the 
coefficients which appear in the short-time expansion of the heat-kernel K(t) correspond- 
ing to a Laplacian-like operator on a D-dimensional manifold M.. In mathematics this 
interest stems, in particular, from the well-known connections that exist between the 
heat-equation and the Atiyah-Singer index theorem In physics, the importance of 
that expansion is notorious in different domains of quantum field theory, where it is com- 
monly known as the (integrated) Schwinger-De Witt proper-time expansion ||, In this 
context, the heat-equation for an elliptic (in general pseudoelliptic) differential operator 
P and the corresponding zeta function (p(s) has been realized to be a particularly useful 
tool for the determination of effective actions |4[] and for the calculation of vacuum or 
Casimir energies || (a fundamental issue for understanding the vacuum structure of a 
quantum field theory). Here usually the derivative Cp(0) °f the zeta function Q and its 
value at s = —1/2 (sometimes the principal part) are needed [|| [J. 

In this paper we would like to exploit another property of the zeta function (p(s) 
corresponding to an elliptic operator P, namely its well-known close connection with 
the heat-kernel expansion. In spite of the fact that almost everybody is aware of such 
connection, its actual use in the literature has remained very scarce till now. If the 
manifold M. has a boundary dA4, the coefficients B n in the short-time expansion have 
both a volume and a boundary part [p], || . It is usual to write this expansion in the form 



K{t) ~ (4vrt)" 



oc 



/ , B k t 

=0,1/2,1,... 



with 



Bi 



dV b n + 



M 



dS 



dM 



;i.2) 



For the volume part very effective systematic schemes have been developed (see for ex- 
ample [[], [TO], |TTJ). The calculation of c n , however, is in general more difficult. Only quite 



recently has the coefficient C2 for Dirichlet and for Neumann boundary conditions been 
found [12-17]. Very new results on the coefficient B5/2 for manifolds with totally geodesic 
boundaries will be given in |18| . 

When using the general formalism of Ref. [12] for higher-spin particles, Moss and 
Poletti [19], |20| found a discrepancy with the direct calculations of D'Eath and Esposito 
21] (see also [22-25]). The latter results have been confirmed in [|6|, 27| , where a new 



systematic scheme for the calculation of C2 has been designed in the context of the Hartle- 
Hawking wave-function of the universe and for the case when the whole set of basis 
functions is known [F2^, [37]. Finally, very recently the discrepancy has been resolved 
completely p8| and now the results that are found using the general algorithm |29] are in 



agreement with those coming from the direct calculations [21-27]. 

The connection between the heat-kernel expansion, Eq. (|1 . 1|) and the associated zeta 
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function is established through the formulas p0 | 



Res COO = 2 (1.3) 
(47r) 2 T{s) 



for s = f , ^i, §; -^±i, for I E 1N , and 



C(-p) = (-W^S, (i.D 

(4tt) ~ 
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for p e INo- The aim of the present article is to show that these equations, ( p..3| ) and ( |1.4| ), 
can actually serve as a very convenient starting point for the calculation of th coefficients 
-8^, even in the cases when the eigenvalues of the operator P under consideration are 
not known. The good knowledge in explicit zeta-function evaluations that have been 
accumulated in the past few years (for a review of many results in this respect, see |]31|) 
will allow us to elaborate a very competitive method of calculation of the heat-kernel 
coefficients which makes use of rather common ingredients, such as integral representations 
of the spectral sum, Mellin transforms, non-trivial commutation of series and integrals 
and skilful analytic continuation of zeta functions on the complex plane. 

To explain the method in detail we will consider the Laplace operator on the D- 
dimensional ball with Dirichlet, Neumann or (in general) Robin boundary conditions. 
Earlier investigations on the first few coefficients are due, for D = 1, to Stewartson and 
Waechter |2|, to Waechter in D = 2 g| and to Kennedy @, |35[ in up to D = 5 



dimensions (for recent results on the functional determinant of the Laplace operator on 



the three and four dimensional ball see [36]). In these references the method was based 



on the use of Laplace transformations of the heat-kernel K{t) itself. In that method an 
intermediate cut off has to be introduced at some point — because one needs to consider 
the Laplace transform of a function which is singular at t — 0. In contrast, in our approach 
it is the complex argument s of the zeta function of the Laplace operator which very neatly 



serves for the regularization of all sums (in just the usual way f31~H ). 

The layout of the paper is as follows. In section 2 we briefly describe the eigenvalue 
problem of the massive Laplace operator on the ball and derive a representation of the 
associated zeta function in terms of a contour integral. We consider the massive Laplace 
operator because the analytical continuation procedure is slightly easier for the case of 
non-vanishing mass. In section 3 we describe how an analytical representation of the 
zeta function — valid in the strip (1 — N)/2 < 9fts < 1 — can be obtained for any N, 
restricting our considerations in this section to D = 3 and to the case of Dirichlet boundary 
conditions. This representation will display very clearly the meromorphic structure of the 
zeta function. As is then shown in section 4, from this representation it is quite immediate 
to read off special properties, as the ones reflected by ( |1.3| ) and ( |1.4| ), in order to find the 



heat-kernel coefficients. In section 5 we explain the small changes in the procedure that 
are necessary in order to treat Robin boundary conditions, in general. Finally, in section 6 
we study the modification to be introduced in the formulas for considering any arbitrary 
dimension D. In appendix A we exhibit some technical details of the calculation and 
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in Apps. B, C and D we give explicit tables of the heat-kernel coefficients for Dirichlet, 
Neumann and general Robin boundary conditions, for the dimensions D = 3,4, 5. 



2 Heat-kernel coefficients on the D- dimensional ball 

As explained in the introduction, we are interested in the zeta function of the operator 
(—A + m 2 ) on the D-dimensional ball B D = {x 6 IR D ; \x\ < R} endowed with Dirichlet, 
Neumann or Robin boundary conditions. The zeta function is formally defined as 

c(«)=£*r> (2.i) 

k 

with the eigenvalues X k being determined through 

(_A + m 2 )0 fc (a;) = \ k <p k {x) (2.2) 

(k is in general a multiindex here), together with one of the three boundary conditions 
above. It is convenient to introduce a spherical coordinate basis, with r = \x\ and D — 1 
angles Q = D _ 2 , V 9 )- m these coordinates, a complete set of solutions of Eq. Q2.2|) 

together with one of the mentioned boundary conditions may be given in the form 

<f>i,m,n(r,Sl) = r 1 '^ J l+ D^2(wi )n r)Y l+ D(n), (2.3) 



with Jj+(£)_2)/2 being Bessel functions and Yi +D / 2 hyperspherical harmonics |j57| . The wi t . 
(> 0) are determined through the boundary conditions by 



J l ! d-2 (wi r,R) = 0, for Dirichlet boundary conditions, 



u 



—J l+ D^(wi n R)+wi n J' l D-2(wi n r)\ r= R = 0, for Robin boundary conditions. (2.4) 

R 2 ' ' ' ~ 2 ' 

As is clear, the case u = (l—D/2) of the (general) Robin boundary conditions corresponds 
to the Neumann boundary conditions. In this notations, using X^ n = wf n + m 2 , the zeta 
function can be given in the form 

CO CO 

c(*) = EE^)Kn+™ 2 r, (2.5) 

n=0 Z=0 

where w^ n (> 0) is defined as the n-th root of the 1-th equation. Here the sum over n is 
extended over all possible roots w^ n on the positive real axis, and di(D) is the number of 
independent harmonic polynomials, which defines the degeneracy of each value of I and 
n in D dimensions. Explicitly, 

d i (D) = (2Z + D-2) ( | ! |^_~^ ! . (2.6) 

Furthermore, here and in what follows the prime will always mean derivative of the 
function with respect to its argument. 
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To distinguish in the notation among the different cases, we will use the indices D, N 
and R to denote Dirichlet, Neumann and Robin boundary conditions, respectively. Thus, 
we will write (d, (n and for the corresponding zeta functions. Using for the moment 
the unified notation $i + (D-2)/2{wi jn R) = for the boundary condition Eq. (|2.4| ), it turns 
out that Eq. ( |2.5| ) may be written under the form of a contour integral on the complex 
plane, 

00 r dk 8 

C( S ) = X>(L>) J — (fc 2 + m 2 )--ln^ +v (A; J R), (2.7) 

where the contour 7 runs counterclockwise and must enclose all the solutions of ( [2. 4]) on 
the positive real axis (for a similar treatment of the zeta function as a contour integral 
, |3l|). This representation of the zeta function in terms of a contour integral 
around some circuit 7 on the complex plane, Eq. (|2.7D, is the first step of our procedure. 

Depending on the value of the dimension D and on the boundary conditions chosen, 
the analysis of the zeta function, Eq. (|2~7] ) — to be given below — will differ, but just 
in small details. For this reason, we will only describe at length the case of the three- 
dimensional ball with Dirichlet boundary condition. The derivation of the analogous 
results for the other boundary conditions and higher dimensions will then be clear, and 
shall be indicated only briefly. 



3 A quick procedure for calculating heat-kernel co- 
efficients 

As explained above, we will illustrate the procedure in the case of the three-dimensional 
ball with Dirichlet boundary conditions. For D = 3 the degeneracy is di(3) = 21 + 1, 
so that the starting point of the calculation reads (we omit further indication of the 
dimension in the notation) 

00 r dk 8 

CdOO = E(2i + 1) J ^ (k 2 + m 2 )- s - In J l+ ,(kR). (3.1) 

i— ■y 

As it stands, the representation ( ^.1[ ) is valid for 3?s > 3/2. However, we are interested in 
the properties of (d(s) in the range < and thus, we need to perform the analytical 
continuation to the left domain of the complex plane. Before considering in detail the 
/-summation, we will first proceed with the /c-integral alone. 

The first specific idea is to shift the integration contour and place it along the imaginary 
axis. In order to avoid contributions coming from the origin k = 0, we will consider (with 
v = I + 1/2) the expression 

Cd = l^n {k2 + ln (*"" UkR) ) ' (3 ' 2) 

7 
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where the additional factor k v in the logarithm does not change the result, for no addi- 
tional pole is enclosed. One then easily obtains 



sm 7TS 



^ fdk [k 2 -m 2 ]- 3 -^- In (k- u UkR) 
it J ok v 



D 



(3.3) 



valid in the strip 1/2 < 3?s < 1. A similar representation valid for m = has been given 
in 



40, 41 



As the second step of our method, we make use of the uniform expansion of the Bessel 



function I u (k) for v — > oo as z = k/u fixed [42]. One has 



Uvz) ~ 



i + E 



U k (t) 



k=l 



(3.4) 



with t = 1 j\J 1 + z 2 and rj = \/l + z 2 + ln[z/(l + yl + z 2 )}. The first few coefficients are 



listed in f42|| , higher coefficients are immediate to obtain by using the recursion |42 

1 1 f 

o 



(3.5) 



starting with Uo(t) = 1. As is clear, all the u k (t) are polynomials in t. Furthermore, the 
coefficients D n {t) defined by 



In 



1 + E 



u k (t) 



k=l 



^ Dn(t) 



(3.6) 



are easily found with the help of a simple computer program. 

Now comes what can be considered as the third step of our method. By adding and 
subtracting N leading terms of the asymptotic expansion, Eq. ( p.6|) , for v — > oo, Eq. (|3~3] ) 
may be split into the following pieces 



N 



Co = Z D {s) + J2 ^> D (s 



with the definitions 

Z»{s) - 



sin its 



71 



dz 



mR/u 



i=-l 



zv\ ■ 



RJ 

-In 



d_ 

dz 



,-,!/(/ 



v / 2^(l + z 2 )i 



| In z v I v (zv) 



Dn{t) 
.,n 



(3.7) 



(3i 



71=1 



and 



Sin 7TS 



7T 



G?2 



mi?/ v 



zv\- 



m 



d_ 

dz 



In (z-"er>) , 



(3.9) 



6 



.4 



u,D 



sin (ns 



71 



dz 



mR/ v 
oo 



zv\- 



A 



sin (ns) 



7T 



dz 



mR/ v 



ZV 

~R, 



m 



m 



ln(l + z 1 



(3.10) 
(3.11) 

The essential idea is conveyed here by the fact that the representation ( |3.7|) has the 
following important properties. First, by considering the asymptotics of the integrand in 
Eq. (|3.8|) for z — > mR/u and z — > oo, it can be seen that the function 



d_ 

dz 

d_ 

dz 



Z D (s) 



oo 1 

S(2Z + l)Z l ^{s) 



1=0 



is analytic on the strip (1 — N)/2 < dts < 1. For this reason, it gives no contribution to 
the residue of Cd( s ) i n that strip. Furthermore, for s = —k, k G 1N , k < —1 + N/2, we 
have Z(s) = and, thus, it also yields no contribution to the values of the zeta function 
at these points. Together with Eqs. ( |1.3j ) and ( |1.4| ), this result means that the heat-kernel 
coefficients are just determined by the terms Af(s) with 



E(2/ + 1)4 + ^' D ( S ) 

1=0 



(3.12) 



As they stand, the A?' n {s) in Eqs. (0), (|3~T0|) and ( pTD are well defined on the strip 
1/2 < 3?s < 1 (at least). And we will now show that the analytic continuation in the 
parameter s to the whole of the complex plane, in terms of known functions, can be 



performed. Keeping in mind that Di(t) is a polynomial in t, all the A^ ,D (s) are in fact 
hypergeometric functions, which is seen by means of the basic relation E3l 



2 F 1 (a,b; c; z) 



T(c) 



dt t b -\i -ty-^ii -tz)- a . 



Let us consider first in detail A 
immediately that 



r(6)r(c-6) Jo 

and the corresponding A^^s), Aq(s). One finds 



m 



-2s 



2^ 
1 



Rm 



r (-i) 



T(s) 



m- 2s 2 F l 



l,s; 1, 





1 


1 




2 : 


2 ; " 






1 


" )1 




— m 


mR) 




4 



mR) / 
1 + 



-2s 



V 

m 

2 



/ v 



\mR 



2s , (3.13) 

— s 

, (3.14) 



where in the last equality we have used that 2-^1(0, s; a; x) = (1 — x)~ s . 

The next step is to consider the summation over /. For A v l x (s) this is best done using 
a Mellin-Barnes type integral representation of the hypergeometric functions 



iFi(a,b; c; z) 



Tie) 



T(a)T(b) 2vri 



dt na+mb+tn-t) 

T(c + t) K ; ' 



(3.15) 
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where the contour is such that the poles of T(a + t)T(b + t)/T(c + t) lie to the left of it 
and the poles of r(— t) to the right ||43|| . After interchanging the summation over / and 
the integration in (|3.15| ), the result will be a Hurwitz zeta function, which is defined as 



Ch(s; v) 



£(*+*)- 

1=0 



$ls > 1. 



(3.16) 



However, as is well known, one has to be very careful with this kind of manipulations, 
what has been realized and explained with great detail in E3, [45], f46|. This point is of 



crucial importance (it has been the source of many errors in the literature over the past 
ten years ||3l|) and can be considered as the fourth step of our original procedure here. 
Applying the method, as described in the mentioned references, to AP x (s), 



E(2i+i: 

1=0 



ni 



-2s 



2^ 



Rm 



r (-§) 



T(s) 



, 1 1 1 

2^i - tt; tt; - 



2' 2' 



mR 



I + i 



it turns out that we may interchange the J2i an d the integral in Eq. (ft.!5|) only if for the 
real part 3?C of the contour the condition < — 1 is satisfied. However, the argument 
r(-l/2 + t)T(s - l/2)/r(l/2 + t) has a pole at t = 1/2. Thus the contour C coming 
from — ioo must cross the real axis to the right of t — 1/2, and then once more between 
and 1/2 (in order that the pole t — of r(— t) lies to the right of it), before going to 
+ioo. That is, before interchanging the sum and the integral we have to shift the contour 
C over the pole at t = 1/2 to the left, cancelling the (potentially divergent) second piece 
in A^^s). Closing then the contour to the left, we end up with the following expression 
in terms of Hurwitz zeta functions 



A^(s) 



d2s oo 
-ft ^ 



-IV o T (j+ S ~l) 



s+j 



Ch(2j' + 2s-2;1/2). (3.17) 



For Aq one only needs to use the binomial expansion in order to find 



£ { —^L(mR) 2 iT(s + j)C H {2j + 2,-1; 1/2). 



2T s 



3=0 



J- 



(3.18) 



The series are convergent for \mR\ < 1/2. These representations (|3.17|) and (|3.18|) show 



very clearly the analytic structure of A^±(s) and A^s). As the fifth (and final) step of 
our procedure, we are left with the quite simple task of explictly evaluating this analytic 
structure, namely of finding its poles and some point values, and of adding all contributions 
together. 

The point values A® 1Q (— p), p G 1N — respectively their residues in s = 1/2, — {21 + 
l)/2, / G INo — necessary for the calculation of the associated heat-kernel coefficients are 



S 



immediate to obtain, using just 

£,(1 + 6,1/2) 

T(e - n) 



1 (-1) 
e n\ 



0( t 



(3.19) 



However, before we can actually calculate (an in principle arbitrary number of) the heat- 
kernel coefficients, we need to obtain analytic expressions for the Af(s), i G IN. As is 
easy to see, they are similar to the ones for A^^s) and Aq(s) above. We need to recall 
only that Di(t), Eq. (|3.6|), is a polynomial in t, 



i+2a 



(3.20) 



a=0 



which coefficients x^ a are easily found by using Eqs. (|3~5| ) and (|37q ) directly, or either 
by using the direct recursion relations presented in appendix A. Thus the calculation of 
A"' D (s) is essentially solved through the identity 



oo 


~(zv\ 2 2 


I ^ z 


— — m 


\rJ 


R/v 





— 

dz 



— m 



n r ( g + f)r(i- s ) 
2{mRy r (i + f ) 



XV" 



1 + 



V \ 



mR 



) 



(3.21; 



The remaining sum may be done as mentioned for A®, and we end up with 

A? {*) = 

1 \?) j=o J'- 

+ 2a)T (s + a + j + 



E ^r(™R) 2j ( H (-l +i + 2j + 2s; 1/2) 



X E/ 

a=0 



r(i 



a + 



(3.22) 



In summary we have obtained the analytic expression of all the asymptotic terms 
coming from expansion (|3.4|) in its most elementary form, which involves the very familiar 
Hurwitz zeta functions and Gamma functions only. Expressions ( |3.17[ ), ( |3.18| ) and ( |3.2z 
constitute the explicit starting point for the calculation of an — in principle arbitrary 
number of heat-kernel coefficients in an extremely quick way. 



4 Heat-kernel coefficients for Dirichlet boundary con- 
ditions on the three-dimensional ball 

Let us now see how the analysis in Sect. 2 can be used for a very effective calculation 
of the heat-kernel coefficients. The dependence of the coefficients on the mass is already 
contained in the coefficients of the massless case through 

K m (t) = K m=0 (t)e- m2t 
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and for this reason we shall restrict ourselves to m = 0. For the sums in ( p.!7|) , ( |3.18| ) 
and ( [3.22j ) this means that only j = will contribute. 

We shall distinguish between the coefficients B& with integer and half-integer index k, 
because the situation is actually different in both cases. In fact, corresponding to Eq. (|1.3|) 
(resp. Eq. ( |1.4j )), the residue of (resp. the value of the function) £c is needed. 

Let us start with the case of integer index k G IN, so that Res Co (3/2 — k) is to 
be calculated. In order that Zjj(s) does not contribute, one has to choose N = 2k — 1 
and thus only the asymptotic terms A?(s), j = — 1, 0, 1, 2k — 1, will provide some 
contribution. Furthermore, one may see very easily which terms in the different A®(s) 
contribute. An important feature is, that for i = 2n, n G INo, Af(s) does not contribute 
to Bk for k G IN. The relevant residues are found to be 

Res A D 1 (--k] = R "T N Or f 1 - 2fc; - 

~ X V2 ) (k -l)!2v/?r (| -k) V '2 

„ An /3 ,\ R 3 ~ 2k ^ (2A;-l + 2a)a! 

Res Ag_ x (■z-k) = -—j- pr £ ^3fc-i.«- 



~ 2 > 2r(|-fc) a = ,a r(l + a + k) ' 

and for n E IN, n <k — 1, k < 3n, 

„ A /3 \ (-l)^ 3 " 2 * / , l\*^ n (-l) a+n (2n + 2a- 1) 

Res A 2n „! [ - - k) = ^-/- ^(h (l + 2n-2k;-) £ x a ^> ± > 

\2 J r(|-ife) V 2/ ^ (fe-i-o-n)! 

whereas for n < k — 1, A; > 3n, we have 

3 ,\ (-l) k R 3 - 2 * / „ n , lX 2 ^ 1 (-l) a+n (2n + 2a- 1) 

n— l,a " 



/3 \ (—1) it / 1\ 

Res A 2n _! fc = - — j- ^-Cu ( 1 + 2n - 2fc; - ] V x 2n - 1 „ 

V2 / r (|_fc) Sii V '2/ a t^ ' (Jfe-1 

From these results we readily obtain the heat-kernel coefficients through 

B k 



n) 



Res Cd (- - fc) = Res (3 " fc 



(47r)fr §-Jfc 



The coefficients up to B w are listed in appendix B. 

Let us now consider the calculation of the coefficients corresponding to half-integer 
index B k+ i/ 2 , k G IN. Here the value of Co (3/2 — k) is needed and one finds N = 2k. It 
is apparent that the Af(s) with odd i, % = 2j — 1, j G IN , do not contribute now. The 
relevant values of the Af(s) read 

T>2-2k , I 

K(l-k) = — CH (l-2k-- 

^ (a + *-!)! 
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and for n £ IN, n < k — 1, k < 3n — 1, 

/ ix fe— n— 1 ^ -| ^n+a 

<(l-fc) = -2i? 2 - 2fe (fc-l)!C^ l + 2n-2£;;- £ x 2n , a - 



2/ (fc-w-a- l)!(a + n- 1)! 

whereas for n < k — 1, k > 3n — 1, we have 

/ 1 \ 2n C— 11 n+a 

<(1 - fc) = -2i? 2 " 2 ^ - 1 + 2. - 2fc; - x 2n , fl 



a=0 



(A;-n-o-l)!(a + n-l)!' 



And from these results, we finally obtain 



k f_l 'jfc— l/t. _ i M 

Cd(1 - fc) = £ ^(1 - k) = 1 J \, J ' i? fc+ i. 
n=o (4tt) 2 



The heat-kernel coefficients are listed in appendix B too. Using (H{—n;q) = 

—B n+ i(q)/{n + 1), n £ INo, the results might have been given, equivalently, in terms 
of Bernoulli polynomials B n+ i(q). 



5 Robin boundary conditions on the three-dimensional 
ball 

When Robin boundary conditions are imposed, using the same method of the preceding 
sections we can write the zeta function as 



C r\k <~l r?/ 

0^) = £(2/ + 1) J — (P + m 2 r^ln [^J l+ i(kR) + kJ' l+ ,(kR) 

i— ty 

and, in analogy with Eq. (|3.3|), we then consider 



d 



(5-1) 



c R = s -^jdk { e- m r s l k in 



u 



k-"{-I v {kR) + 



(5.2) 



Employing the same idea as for Dirichlet boundary conditions, this time we have in 
addition the following uniform asymptotic expansion 



I'(uz) ~ 



1 e^(l + z 2 U 



\j2ixv z 
with the Vk{t) determined by 

Vkif) =u k (t)+t(t 2 -1) 



1 + f 5W 



k=l 



(5.3) 
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In analogy with Eq. (|3.6|), we write 
In 



fc=i ^ v v S ^ 



n=l 



where the functions M n (t) are easily obtained. At this point we see already, that for Robin 
boundary conditions no additional calculation is necessary. Comparing the expansion 
( |5.3| ) with ( |3.4| ) and introducing Af(s) for the contributions coming from the asymptotic 
terms, one has 

A* 1 (s)=A D 1 (s), A*( S )=-A°(s). (5.5) 
Furthermore, the functions M,(t) are of the form 

M t (t) = jrzi,at l+a (5.6) 

a=0 

(notice that here, in contrast with the case of Dirichlet boundary conditions, all powers 
between i and 3i are present). As a result, we find 

1 V s ) j=o 3- 

v ( i + a ) r ( s + j + i±2) 
x> z ia ^ r '-. (5.7) 

One can show again that only the even indices i contribute to the residues of (r(s), 
whereas the odd ones will contribute to the point values. 

Restricting ourselves as before (see the comment in the previous section) to the mass- 
less case, the results for the heat-kernel coefficients may now be read off from the formulas 
in the previous section. One has 

Res A\ (- - k^j = Res A° x (J - kj , 



/ q \ / j "iK— l.o /i 

2r(|-*)S r(| + /, + f 



where the expressions for A 2n _ x are found from the results in Sect. 3, once x 2 n-i,a has 
been replaced with Z2n-i,2a- The coefficients for Neumann boundary conditions are given 
in appendix C, and for the general case (u arbitrary) in appendix D. For the point values 
the analogous formulas read 

A«(l-k) = -Ag(l-k), 

4fc (2k + a)T ( 1 + %) 
A? k (l-k) = (-l)^ 2 ->-l)!£W ^ 2) 



a=0 



2T (l + k + f 



and once more the replacement of X2 n ,a with z 2 n,2a leads to the results for A 2n - The results 
for the heat-kernel coefficients are summarized in Apps. C and D. 
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6 Generalization to the D- dimensional ball 



As we will now explain, for the generalization of our results to the case of a D-dimensional 
ball almost no additional calculations are necessary. Let us discuss first the case of 
Dirichlet boundary condition. The starting point of the analysis is now 







^ r f If 

Cd(s) =Y^d l (D)J — (k 2 + ///-') In J, , n , 



1=0 



dk 



(kR). 



(6. 



It is easy to see that the above treatment for the individual terms of the /-series, 



(6.9) 



remains valid, once we have set v = I + (D — 2)/2. In order to use our procedure for the 
whole /-summation, what remains to be done is to substitute for the degeneracy di(D) 
its value in powers of / + (D — 2)/2, in order to find again expressions in terms of the 
Hurwitz zeta function (h(s; (D — 2)/2). Writing 



D ~ 2 ( D — 2 

di(D)= E e "P) / + — — 

a=l V / 



(6.10) 



the final results for A^^s), Aq(s) and Af(s), i G IN, may be read off from Eqs. (|3.17|) , 
( CT ) and |T22|). We find 



[D-2 



S+j 



(6.11) 



A?(s) 



J2 e a {D)( H {2j +2s-l-a;{D- 2)/2) 

a=l 

_gL_i (mfl) « r(3 + j) 

'D-2 

]T e a (D)( H (2j + 2s - a; (D - 2)/2) 

.a=l 

R 2s oo /_ 1 y 



(6.12) 



2T s 



E 



(mR) 



2.7 



j=0 J- 

D-2 

e a (D)( H (-a + i + 2j + 2s; (D - 2)/2) 



(6.13) 



,«=i 



+ 2a)T(s + a + j + 



r (l + a + f 



We shall spare the reader the analogous results for Robin boundary conditions. They 
need not be given explicitly, since the procedure is absolutely clear by now. Let us just 
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write down the relevant residues and point values of (d(s) (the Robin case follows from 
the replacements explained in Sect. 5). They read 

/ o x p3-2fc D-2 , n _ 9\ 

«(H - ( ;-i)u^(i-.) 5 e ° (gK "( 2 - 2t - a: — )■ 



for n — 1, k — 1, k > 3n, 

R-es A 2n -i - k) = (-l) k f ^ Y^e a (D)( H (2 + 2n-a-2k;^ 2 
/ 2T - kj a =i v 1 

lV , +a (2n + 2a-l) 

X 2^ (-1) X 2n -l,a " 7—^ 

a =o ' (k - 1 - n - a)\F (J + a + 

whereas for < n, it reads 

(o \ p3— 2fc £?— 2 , j~\ 2\ 

2" fc ) = (-l)" 2r ^ _ ^ E e*{D) CH[2 + 2n-a- 2k; —^-j 



^ \ ,xn+a 2n + 2a-l 

X E (- 1 ) ^2n-l,a- " ~ 

(k - 1 - n - a)\T (§ + a + 



For higher indices it is adviceable to distinguish between .D even and D odd. For D odd 
contributions arise for n = k, + (D — 3)/2, and read 

D , /3 \ i? 3 - 2fe 2 ^ (2n + 2a - l)(a + n - k)\ 

Res A 2n -i - — k J — — —77 -ye 1+2 „_2fc 2^ x 2ri-i,a — r: \ , 

4r (I - A;) tl> r (i + a + n) 



whereas for D even the indices run from n = k, k + (D — 4)/2, and the results are 
. / 3 \ i? 3 - 2fe ^ r(|-A; + a + n) 

ReS A 2ri - — k) — j- -Y e 2+2n-2fc y. %2n,a 7 ; TTj • 

V2 J 2r (I - A;) (a + n-1)! 

Let us conclude with the list of point values. The leading asymptotics A D ± gives only 
contributions for k — 0, 



A-i(0) = -\ ^e a (D)C H (-« - 1; ^ 
Furthermore, for n — 1, k — 1, we have 

A 2n (l-A;) = -i? 2 - 2fc (fc-l)! X; 2 ea ( J D)C„(-a + 2n + 2-2A;;^^) 

a=l V 2 ' 

2n 

X a ^ X2 "' a (a + n-l)!(A;-l-a-n)!' 
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if k > 3n — 1, and if k < 3n — 1 



D-2 



A 2 „(l-A;) = -i? 2 ~ 2fc (A;-l)! E e Q p)Cff -« + 2n + 2 - 2A;; 



a=l 



D-2 



k—n—l 



-1) 



a+n 



X £ X2n ' a (a + n-l)!(A;-l-a-r2)!- 



Finally, for Z) odd and for n = k, k + (D — 3)/2 

1 



^,,(1 - />') ^(-l) fc (fc-l)!i? 2 - 2fc e 1+2 „_ 2fc i:x 2n , a -f a '" _/ '' ,! 

a=0 



±^ (a + n - 1)1 

whereas for D even the indices run from n = k + 1, k + (D — 2)/2, and the result reads 



1 



2n-l 



A 2n _! (1 - fc) = -(-1)"(A; - l)\R 2 - 2k e 2n - 2k £ x 2n _ lja 



(2n-l + 2a)r (i - fc + a + n) 



a=0 



r ( | + a + n 



The formulas above simplify a bit if we write the degeneracy (|2.6D under the form 



di{D) 
di(D) 



(D-2)! 
2 



1 + 



D-2 



1 + 



D-2 



D 



£-2 
2 



x ... x Z + 



X ... X / + 



D-2 



D-2 



for D odd, 



for D even, 



so that e 2 fc(D) = for D odd and e 2 fc_i(D) = for D even, G IN. 

Furthermore, one might use the following recursion for the coefficients e a (D) appearing 
in the expression of the degeneracy di(D), Eq. ( |6.1U| ), 

1 



e 2Q (D + 2) 
e 2a _x(D + 2) 



D(D-l) 
1 

D(D-1) 



D \ 2 
e 2a - 2 (D) - ( - - lj e 2Q (D) 



e 2a -3p) - ( y - l) e 2a ^{D) 



for D even, 
, for D odd, 



where we have used the definitions e-k{D) = for k G INq and e a (D) = for a > D — 2. 
We have performed explicit calculations for D = 4 and D = 5. One has in these cases 



d,(4) = (Z + l) 2 , e 1 (4) = 0, e 2 (4) = l, 



d*(5) 



1 / 3 
- Z + - 
3 V 2 



Z + 



ei(5) = -^, e 2 (5)=0, e 3 (5) = I 



The results for the heat-kernel coefficients are presented in Apps. B, C and D. 
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7 Conclusions 



As promised in the introduction, we have developed in this paper a very convenient 
method in order to deal with the problem of the calculation of heat-kernel coefficients 
corresponding to an arbitrary elliptic operator with any of the usual boundary conditions 
(Dirichlet, Neumann or Robin), with the only proviso that the behavior of some basis for 
its spectrum should be known (even if the eigenvalues themselves are actually unknown). 

This is indeed a very common case in mathematical physics, what conferes to our 
procedure a wide generality of application. Another fundamental characteristic of the 
method is its extreme simplicity, which comes in part from the quite strong background 
on zeta function computations that we have acquired during the last half a dozen years. 
This knowledge conferes to the new method the same elegance that the procedure of zeta 
function regularization (including the analytic continuation techniques and non-trivial 
series commutation that it involves) has in itself. 

Finally, we have tried our method with explicit examples and give several tables of 
heat-kernel coefficients that have been calculated here (with relative easiness) for the first 
time. For the near future we envisage to investigate other physical applications where the 
method can prove useful. 



Note: At the final stage of our analysis, P. Gilkey made us aware of related research by 
M. Levitin ||7], who has further developed the approach of Kennedy [|4], [35[, also with the 



aim of calculating higher-order heat-kernel coefficients. We are indebted with M. Levitin 
for sending us his results, which have served as a very good check of our calculations. All 
results in common with his are in complete agreement. 
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A Recursion relation for the coefficients x 



i,a 



In this appendix we present the recursion relations for the coefficients Xi >a , Eq. (|3.20|) . For 
convenience let us introduce for i E IN, a — 0, i, 



c i+l,a 



2 i+1 (i + 2a) 

Then, starting with c^o = — 1, we find the following recursion relation, 

1 

Q,0 — (i ~ 2)02-1,0 — - X] C i-s,0 c s,0, 
^ 8=1 

1 <-1 
1 s=l 

and for a = 1 2 — 2, we have 



Ci,a — {i — 2 + 2a)(cj_i j(I — Cj_i j0 _i) 

^ i-l / Mm(a,s— 1) Afm(a— l,s— 1) 

— 9 I C i~s,a-jC s ,j — ^ Cj_ Sja __j_iC SJ 

s=l \j=Max(0,l+a+s-i) j=Afax(0,a+s-j) 

This relation can be used very effectively for the calculation of the coefficients Xi a . 



B Heat-kernel coefficients for Dirichlet boundary con- 
ditions 

In this appendix we list our results for the heat-kernel coefficients of the Laplace operator 
in 3, 4 and 5 dimensions with Dirichlet boundary conditions. Here and in the following 
appendices, the first coefficients B , ...,B 5 / 2 are listed for completeness and may also be 
found in or derived from [E?|. 



In three dimensions we have found that 

Bo = -tiR 3 

B 1/2 = -2^ 2 R 2 
8irR 



3 



B^ = -U i/2 



'3/2 

B 2 



B5/2 



6 
16tt 

315R 

7T 3 / 2 

120i? 2 



17 



64vT 

3 9009i? 3 
47tt 3 / 2 
7/2 - 2016CLR 4 



B A 



202816tt 
72747675^ 
521tt 3 / 2 



59/2 443520.R 6 



B 5 = 

B 

B 6 = 

B\z/2 = 

B 7 = 

B 

B 8 = 

B 

B 9 = 

B 

Bio = 

In four dimensions the result is 



25426048tt 

15058768725/? 7 

9521tt 3 / 2 

112 ~ ~11531520i? 8 

90878576896tt 



67689165418875i? 9 
34344493^3/2 

~4704860160(LR 10 
22835854180352tt 

'17531493843488625/? 11 
36201091vr 3 / 2 

15/2 ~ ~47048601600i? 12 

1509389910845640704tt 

1019964780320324713875i? 13 
153984929039tt 3 / 2 

17/2 ~ ~ 16448191 119360CLR 14 

1673450232605639069696tt 

872477873086005760248675.R 15 
13334525091737tt 3 / 2 

" 2 ~ 10362360405196800i? 16 

643985013732181345325056tt 

231206636367791526465898875i? 17 ' 



Bo = \* 2 R A 

B 1/2 = 
B 1 = 2n 2 R 2 

llir^R 

B3/2 = - 



B; 

B 



32 
4tt 2 



45 

35tt 5 / 2 
' 2 ~ 4096i? 



18 



464tt 2 
'45045i? 2 

9117T 5 / 2 

196608/? 3 
107456tt 2 

' 14549535/? 4 
827315tt 5 / 2 

201326592i? 5 
23288576tt 2 

'3011753745i? 6 
1585902737T 5 / 2 

11 2 322 122o 172( )/,'• 

200645457927T 2 



B 3 = 

B7/2 = 
B 4 = 

Bg/2 = 

B 5 = 

B 

B, 

Biz/2 

B 7 

B 

B 

B 

B 9 = 

B 

Bio 



1933976154825i? 8 
630648945109tt 5 / 2 

'86586540687360i? 9 
492912963584tt 2 



29464695535275i? 10 
70309732006867tt 5 / 2 

152 ~ ~5541538603991040i? n 
37648078688043008tt 2 



1204208713483264125/? 12 
1578924180477650401tt 5 / 2 

17 2 ~ ~ 62419890835355074560.R 13 
8875043738202275847T 2 

13409327181833639595i? 14 
1018264365864160946171tt 5 / 2 

102 ~ 17976928560582261473280/? 15 
252629551155828479492096tt 2 



1616829624949591094167125.R 16 ' 
Finally, in five dimensions we obtain 

87T 2 i? 5 



B 

B1/2 



15 

47T 5 / 2 ^ 4 

3 

32tt 2 R 3 

Bl = 
B 3/2 = -r >/2 R 2 
128% 2 R 
945 



Bo = 



19 



B 5/2 



B 3 

B 

B 4 
B, 

B 5 
Bn/2 
B 6 



B13/2 



B 7 

B 

B 

B 

B 9 

B 

Bio = 



17tt 5 / 2 
360 

12167T 2 



45045.R 
157tt 5 / 2 

7/2 30240i? 2 
235264tt 2 

43648605/? 3 

5tt 5 / 2 

" 2 21(i l/, , ; 

779264tt 2 



280598175i? 5 
593tt 5 / 2 

449280i? 6 
91757946368tt 2 

43074923448375.R 7 
32815499tt 5 / 2 

28229160960i? 8 
22103738934272tt 2 

105188963060931757? 9 
1 190343 19tt 5 / 2 

' 71 ~ 94097203200/? 10 

533003666 10079744tt 2 

21397862524202616375.R 11 
7986089796017T 5 / 2 

17/2 " 493445733580800i? 12 

381809787573414866944tt 2 

111856137575128943621625/? 13 
146801666871373tt 5 / 2 

19/2 ~ 62174162431180800i? 14 

318152827895794391 12192tt 2 

6031477470464126777371275i? 15 



C Heat-kernel coefficients for Neumann boundary con- 
ditions 

Here is a list of the results we have obtained for the heat-kernel coefficients of the Laplace 
operator in 3, 4 and 5 dimensions with Neumann boundary conditions. In three dimensions 
we have found 



Bq = -;t/, ,:! 



3 



20 



In four dimensions 



B 1/2 = 2n 3 / 2 R 2 
SttR 

Bl ~ ~T 

B3/2 = gvr 3/2 

167T 

~9R 

47tt 3 / 2 

60R 2 
6464tt 



B 



2 



-65/2 



B> 
B 4 



6435i? 3 
3973tt 3 / 2 

7 1 " 10080.R 4 

14766656tt 



5 5 

Z3 
Z? 

fig 

B 

Bio 



31177575.R 5 
5057vr 3 / 2 



59/2 2816CLR 6 



2314167424tt 

10756263375i? 7 

2320069tt 3 / 2 

1 ; 2 ~ 27675648~R8 

1439468204288tt 



13537833083775i? 9 
112984728317T 3 / 2 

13/2 ~ 250925875200i? 10 

369968178163712tt 

5843831281162875i? n 
1718717967893tt 3 / 2 

15/2 57211099545600i? 12 

48366532825354366976tt 

1019964780320324713875.R 13 
113384991528329tt 3 / 2 

17/2 451150384988160I1R 14 

781980237125923045376tt 

17805670879306240005075i? 15 
33839928581307889tt 3 / 2 

11,2 1326382131865190400i? 16 

143924362 16775440050663424tt 

297265675330017676884727125i? 17 



Bo = ^rr 2 R 4 



1 

2' 



21 



B 1/2 = ^' 2 R 
B x = 2n 2 R 2 



3 



B3/2 = 

B 2 = 

B, 

B. 



32 

1167T 2 



45 
5861tt 5 / 2 

'"' 2 ~ 4096i? 

99472tt 2 



B 5 = 

B 

B 

B 

B 

B 

B 

B 

B. 

B 

B\o = 

And, finally, in five dimensions 

B n = 



45045i? 2 
388657vr 5 / 2 

7/2 393216i? 3 
18334144tt 2 

14549535^ 
_ 91095533tt 5 / 2 

9/2 " 201326592i? 5 
6269294336tt 



B A = 



2 



15058768725/? 6 
2096614963tt 5 / 2 

lV2 32212254720/? 7 

1448614636544tt 2 

13537833083775i? 8 
130411491766317T 5 / 2 

1:1/2 ~ ~ 86586540687360i? 9 
38509398708224tt 2 

100179964819935/? 10 
14987877600614637T 5 / 2 

15/2 ~ 5541538603991040i? n 

7562397933317668864tt 2 

13246295848315905375/? 12 
23865356170241004641tt 5 / 2 

17/2 ~ 62419890835355074560i? 13 
3004505191361 1575296tt 2 

36622112051226326625.R 14 
135252966433194092697787tt 5 / 2 

m/2 ~ ~ 233700071287569399152640.R 15 
307843753219621367054336tt 2 

' 230975660707084442023875i? 16 



8tt 2 R 5 
15 



22 



3 

32ir 2 R* 
9 

3n 5 / 2 R 2 
1024ir 2 R 

135 
1873tt 5 / 2 

360 
63296tt 2 

6435i? 

101217T 5 / 2 

1890/? 2 
5040647T 2 

61047i? 3 
198463tt 5 / 2 

55440i? 4 
125689856tt 2 

30879225/? 5 
34154807tt 5 / 2 

34594560i? 6 
564471 70574848tt 2 

~ 157941385977375i? 7 
16602940093vr 5 / 2 

14114580480/? 8 
945576485184512tt 2 

~281253911927625 J R 9 
13550828636809tt 5/2 

5721109954560/? 10 
25910401 1527854628864tt 2 

~ 55634442562926802575i? n 
53795807052692597T 5 / 2 

~ 1973782934323200i? 12 
46180677500935662030848tt 2 

~9587668935011052310425i? 13 
2640354677256557617tt 5 / 2 

994786598898892800/? 14 
1401638457879249954799616tt 2 

" 3067757227347963641 74125i? 15 
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D Heat-kernel coefficients for Robin boundary con- 
ditions 



We conclude our list of results with the leading coefficients for general Robin boundary 
conditions for D = 3, 4 and 5. 

In three dimensions, we have found 

Bo = ^ 
3 

B 1/2 = 2^ 2 R 2 

3/2 

B 3/2 = ^(1 + 24m 2 ) 
2tt 

B 2 = (1 - 18m + 60m 2 - 120m 3 ) 

45iT ; 

71-3/2 

B 5/2 = 77^(2- 15m + 60m 2 -120m 3 + 120m 4 ) 
60ii z 



TT 



2 



B s = 4rn<rw 1633 - 12870m + 46904m 
45045/? 3 

-107536m 3 + 144144m 4 - 96096m 5 ) 

71-3/2 

^ = I»< 151 - 1008,1 + 3612 " 2 

-8400m 3 + 13440m 4 - 13440m 5 + 6720m 6 ) 

B A = -(8243319 - 51363270m + 169826940m 2 

436486050i? 5 V 

-395830040m 3 + 676878800m 4 - 835097120m 5 
+665121600m 6 - 266048640m 7 ) 

TT 3 / 2 

59/2 = 1774080^ ( 14639 ~ 80784m + 249304m 
-556600m 3 + 976800m 4 - 1330560m 5 
+ 1340416m 6 - 887040m 7 + 295680m 8 ) 

B 5 = -(3517532467- 17760354570m 

301175374500i? 7V 

+49945523040M 2 - 105573378240m 3 + 182023225440m 4 
-259648898880m 5 + 295543449600m 6 - 252181862400m 7 
+142779436800m 8 - 40794124800m 9 ) 

In four dimensions, the results read 

R 77 K 

-DO = „ 
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B 1/2 = vr 5 / 2 /? 3 
Bi = -2ir 2 R 2 (l + 2u) 

f 5/2p 

B 3/2 = -33- (9 + 32m + 64m 2 ) 
4vr 2 

£ 2 = - — (1 + 30m 3 ) 
45 



Tj-5/2 

5 5/2 = (59 - 224m + 2048m 3 - 4096m 4 ) 

1 4096iT ; 

B 3 = ^r(75 - 286m + 286m 2 + 858m 3 

45045.R 2 V 

-3003m 4 + 3003m 5 ) 

Tj-5/2 

B 7/2 = -(5807- 21024m + 29952m 2 

7/ 393216i? 3V 

+7168m 3 - 110592m 4 + 196608m 5 - 131072m 6 ) 

327T 2 

B ' = -14549535ll( 11726 - 39368 " + 62016 " 2 
-36176m 3 - 75582m 4 + 230945m 5 

-277134m 6 + 138567m 7 ) 

71-5/2 

B 9/2 = -(2961171 - 9105152m + 14440448m 2 

1 201326592i? 5V 

-13142016m 3 - 458752m 4 + 25427968m 5 
-46137344m 6 + 41943040m 7 - 16777216m 8 ) 

647T 2 

B > = - 15058768725* (M19236 ~ 179766 °°* + 274482 °°" 
-28336920m 3 + 14866740m 4 + 14709420m 5 

-49365705m 6 + 65189475m 7 - 47805615m 8 

+ 15935205m 9 ) 

Finally, in five dimensions we have found 

8n 2 R 5 

Bn = 

15 

47T 5 / 2 i? 4 



Bl/2 = 



3 

87T 2 i? 3 



B 1 = _( 5 + 6m) 

7T 5 / 2 R 2 

B3/2 = -^(3 + 8m + 8m 2 ) 

477-2 p 

B 2 = ——(-5 + 6m -60m 2 -120m 3 ) 



25 



Vl-5/2 

B 5/2 = — (-17 - 240m 2 + 480m 4 ) 

27T 2 

= I5Fl35S( 87+ 13442 "- 35464 " 2 



+61776m 3 + 48048m 4 - 96096m 5 ) 



TT 



5/2 



B 7/2 = „ 9 (-88 + 483m - 1806m 2 + 2940m 3 

7560/t 2 

-1680m 4 - 3360m 5 + 3360m 6 ) 

7T 2 

B A = -(-539501 + 4050078m- 12086660m 2 

43648605/? 3 V 

+23878744m 3 - 23715952m 4 + 1478048m 5 
+26604864m 6 - 17736576m 7 ) 

71-5/2 

B g/2 = —(-18927 + 99616m - 302720m 2 

1 2661120 J R 4V 

+576048m 3 - 748704m 4 + 473088m 5 
+ 177408m 6 - 591360m 7 + 295680m 8 ) 

Bs = 90352 6 ' r i2350^ ( - 935536567 + 4964319990 " 
-13111462800m 2 + 25019918880m 3 

-34365190560m 4 + 32451298368m 5 

-12409401600m 6 - 12609093120m 7 

+20397062400m 8 - 8158824960m 9 ). 

This concludes our lists of explicit tables for the heat-kernel coefficients. In the same way, 
results for any desired dimension D are very easy to obtain from the formulas in the text. 
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